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ABSTRACT 

We propose a proper method to measure the size of the narrow hne region 
(NLR) in distant quasars. The apparent angular size of the NLR is, in general, 
too small to resolve technically. However, it is possible to map the NLR if with 
gravitational lensing. In our method, we directly compare the observed image of 
the NLR with the expected lensed images of the NLR for various source sizes and 
lens models. Seeking the best fit image via the comparison procedures, we can 
obtain the best-fit size and the best-fit lens model. We apply this method to the 
two-dimensional spectroscopic data of a famous lensed quasar, Q2237-I-0305. If 
the lens galaxy resembles the applied lens model, an upper limit to the NLR size 
can be set 750 pc. Further, we examine how the fitting results will be improved by 
future observations, taking into account the realistic observational effects, such 
as seeing. Future observations will provide us more stringent constraints on the 
size of the NLR and on the density profile of the lens galaxy. 

Subject headings: galaxies: quasars: emission lines, galaxies: structure, gravita- 
tional lensing 



1. Introduction 

It is widely believed that Seyfert galaxies and quasars have a rather complex structure 
in their nuclei. The central engine seems to be a combination of a supermassive black hole 
and an accretion disk, which are surrounded by the broad- line regions (BLRs), dusty tori. 
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and the narrow line regions (NLRs) (e.g., Antonucci 1993). Therefore, investigating the 
nature of the NLR is useful to obtain good insight into the physics of the central engine and 
its environments. So far, various kinds of observational programs have been carried out to 
understand the basic physical processes involved with the gas in the NLRs, such as photo- 
ionization, shock excitation, and kinematics. Schmitt et al. (2003), for example, performed 
narrow band imaging observations of near-by Seyfert galaxies by using spatially extended 
[O III] A 5007 emission line and measured the size of the NLRs. They found a correlation 
between the size (-Rnlr) and the luminosity (L[o mj) of the NLRs as -Rnlr oc L?q '^j^j? °^, 



which significantly differs from that of quasars (Bennert et al. 2002); -Rnlr oc L^'^^^^j 
On the other hand, Netzer et al. (2004) performed slit spectroscopy of distant quasars, and 
claimed the existence of two distinct populations for luminous active galactic nuclei; that is, 
quasars with and without the NLR. We here note, however, that such arguments are purely 
based on the empirical relation between i?NLR and L\p m] obtained from a relatively small 
number of samples (Bennert ct al. 2002). It is not yet clear if this relation can be accurately 
extrapolated to luminous active galactic nuclei and/or high-redshift quasars. Therefore, 
direct measurements of the size of the NLRs in luminous and/or high-redshift quasars are 
necessary to confirm their findings and to explore the underlying physics. However, it is not 
easy to make direct measurements for distance sources, since the more distant the quasars 
are, the smaller the apparent size of the NLRs becomes. To make matters worse, it is hard 
to measure the size from narrow band imaging observations, because the wavelength of the 
emission line is redshifted and move out from the passband of existing narrow band filters. 

We, here, consider an alternative method which utilizes the gravitational lensing, since 
the gravitational lensing can spatially stretch the source image. In principle, the size of 
the NLRs in multiple quasars can be measured more precisely, compared with that of the 
unlensed quasars with the same intrinsic luminosity. Thus, there have already been such 
attempts by using a two-dimensional spectrograph attached to a ground based telescope 
(Adam ct al. 1989). 

Mcdiavilla et al. (1998) performed two-dimensional spectroscopic observations of a 
lensed quasar with quadruple image, Q2237+0305 (or Huchra's lens), by using an opti- 
cal fiber spectrograph, INTEGRAL. They focused on an arc-like feature in the image of not 
[O III] A 5007 but C III] A 1909 emission hne and derived the size of the NLR in Q2237+0305 
to be ~ 400/i~^ pc. Motta, Mediavilla, Munoz, Falco (2004) used the same observational 
data, but decomposed the spectra into two components; the broad line component and the 
narrow line component. They then made two maps for the two spectral components, finding 
a significant difference between them in terms of the spatial extent of the line emitting re- 
gions. After making two images, they focused on the arc-like feature again and derive a size 
of the NLR to be 700 ~ 900 pc. This value is consistent with that obtained by Mediavilla et 
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al. (1998). Wayth, O'Dowd, Webster (2005), in contrast, could not find an arc-like feature in 
the data taken by the GMOS Integral Field Unit (IFU), claiming that the size of the NLRs is 
significantly smaller than that obtained by Mediavilla et al. (1998). They remove the seeing 
effects via the deconvolution technique, but did not decompose the spectra into the narrow 
line and the broad line component. Since their analysis are different not only in the used 
data but also in the adopted analysis procedures, we cannot simply compare their results 
nor conclude which is more appropriate. Importantly, the lens model degeneracy, which was 
pointed out by Wambsganss, Paczynski (1994), was not properly considered in these works. 

In this paper, we consider a more reliable method to measure the size of the NLRs in 
lenscd quasars, taking into account the effect of seeing and the lens model degeneracy in a 
proper way and then apply our method to both real data and simulated data. In the next 
section, we explain our method to measure the size of the NLRs. The results obtained by 
the currently available data is presented in section 3. In section 4, we examine a potentiality 
of our method for future observational data. The final section is devoted for concluding 
remarks. A concordance cosmology with flm — 0.3, fl\ — 0.7, and Hq — 70 km s~^ Mpc~^ 
(or h — 0.7) is used throughout this paper. 



2. Our Method 
2.1. Procedures 

For a given emissivity profile of the source and for a given lens model, we can calculate an 
ideal image of the extended lenscd source (sec next subsection for the lens model). Convolving 
the ideal image with a point spread function which represents the seeing effect, we can obtain 
an expected image for the observation with an infinite spatial resolution. Re-sampling this 
expected image with a finite spatial resolution, which is same as the spatial sampling rate of 
an observational instrument, we can produce an image expected by the given observational 
instrument. By calculating such images for various source sizes and various lens models, and 
comparing them with the observed one, we can finally obtain the best-fit source size and the 
best-fit lens model. 

Here, we apply x^— minimization method to find the best-fit solution. The total 
is calculated as, 

^ohs / p p \ 2 

2 \Johs,i /model,v /-i \ 

X - -2 > 

i=i » 

where A^obs, /obs,?, /modei,i, ^ud CTj are the number of data points, the observed fiux of the 
i-th data point, the i-th predicted fiux from a given model, and the observational error of 
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the i-th data point, respectively. 

Unless otherwise specified, we assume that the source is circularly symmetric and that 
the emissivity profile of the source (e) is expressed by Gaussian; 

e(r) oc exp (—^r-) ' (2) 

where r represents the distance from the source center. 

Further, we assume that the point spread function due to the seeing effect is also cir- 
cularly symmetric and has a Gaussian profile with the same "full width half maximum" as 
that of the seeing size. 



2.2. Lens model 

It is too hard to measure the size of the NLR without any constraint on the lens model, 
we put a practical assumption on the model that the mass density profile of the lens galaxy 
is elliptical and is expressed by a power-law in radial profile. Since there exists the so-called 
lens model degeneracy (e.g., Wambsganss, Paczynski 1994) even in such situation, we need 
lens models with various density profiles to find reliable constraints on the size of the NLR. 
For simplicity, we adopt so-called "generalized pseudo-isothermal elliptic potential" or "tilted 
Plummer family of elliptic potential" (Blandford, Kochanek 1987). In this model, the lens 
potential (0) at an image position of 6* = {9x, Oy) is given by 

4>0) = 4-^^ K + (1 - e) + (1 + e) el] ' , (3) 

where a^, e, A and 6^ represents the lens size, the ellipticity, the power index and the core 
radius of the lens potential, respectively. There are other models that represent ellipticity 
of the lens galaxies, such as softened power-law eUiptical mass distribution model (Barkana 
1998), or other similar lens models with power- law density profile (Schramm 1990; Keeton 
2001). Even more sophisticated lens models that consist from bulge, disk, and halo compo- 
nent can be used in this kind of studies. It is also interesting to perform the same study by 
using different type of lens models, but main purpose of the current study is to investigate 
how the method works. Moreover, such models require relatively much computational ex- 
pense, and in a current paper, we only focus on a lens model which represented by equation 3. 
Dependence on different type of lens models will be discussed in a future paper. 

A model with large (small) A corresponds to a shallow (steep) density profile. In par- 
ticular, the density profile with A — 0.5 is identical to that of the singular isothermal sphere 
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model. Since the fifth image has not been detected yet (Falco et al. 1996), we set = 
which is equivalent with the lack of the fifth image for this model. As was claimed by Kassi- 
ola, Kovner (1993), e > should be satisfied in this model, since otherwise the iso-density 
contour shows concave shape and such model must be unrealistic for a single lens object. 

For various A, we try to reproduce the positions of the quadruple image of Q2237+0305 
relative to the position of the lens galaxy ^. The number of observable is 2 x 4 = 8, i.e., x— 
and y— coordinate of the four images. Since observed fluxes may be affected by microlensing 
and/or differential dust extinction, we do not take into account the flux ratios between the 
four images in our fitting procedure. The number of free parameters is 6 that is, a^, e, A, x— 
and y— coordinate of the source, and the direction of the major axis of the lens model. Thus, 
degree of freedom in this fitting procedure is 8 — 6 = 2. values of the best-fit lens model 
for various A are presented in the bottom panel of figure 1. The resultant values of cke and 
e which obtained by this fitting procedure is also presented in the middle and the top panel 
of figure 1, respectively. Within 3-a confidence level, a wide range of A, A = 0.41 ~ 0.77, 
is acceptable for the model to reproduce the observed image positions. Even if we restrict 
l-cr confidence level, the model with A = 0.44 ~ 0.73 provides an acceptable fit, and the 
acceptable range does not shrink so much. This represents the lens model degeneracy as was 
found through the previous studies. 

From this fitting procedures, we obtain the best-fit parameters of the lens model for a 
given A as shown in the middle and the top panel of figure 1. Hereafter, we treat the density 
profile (or A) as a single free parameter for the lens model in our proposed method; once 
we change A value, other parameters such as e are automatically replaced with the best-fit 
parameters for new A value. We put a constraint on e value to keep a condition that is 
claimed by Kassiola, Kovner (1993) during the fitting procedures, i.e., e < We can see 
the effect of this constraint in the top panel of figure 1. Basically, the best fit value of cke 
and e becomes lager, when the applied A value for the fitting procedure becomes smaller. 
Though the best fit value of gradually increases with decreasing A until the smallest value 
of A for the fitting (A = 0.4), that of e faces the limiting value of the constraint on the fitting, 
i.e., e = at A ~ 0.46. Below this A value, the best fit value of e is almost identical to 
the limiting value of e. Therefore, we would like to mention that the validity of the fitting 
results at this small A regime is somewhat uncertain. 

Finally, -Rnlr and A is the parameter to be determined from our method. Calculated 
ideal images, which are to be compared with the observed image of Q2237+0305, are shown 
in figure 2. In the case of large i?NLR and/or large A, the expected image tends to exhibit 
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an arc- or ring-like morphology, as is clearly seen in figure 2. Even if we consider the seeing 
effect, obviously, we can recognize differences between the images for the different parameter 
sets. Therefore, we can, in principle, constraint the size of the NLRs and the lens model via 
our proposed method. 

3. Analyses using the Current Data 

Motta, Mediavilla, Munoz, Falco (2004) have performed Gaussian fit to the emission 
line profiles in observed spectra, and have decomposed the detected emission line into a 
broad line component and a narrow line component. The resultant line intensity peak (/) 
and line width (w) of both components of the C III] A 1909 emission line are shown in 
table 1 of Motta, Mediavilla, Munoz, Falco (2004). Since total flux of an emission line with 
Gaussian profile is calculate by y/nlw, we evaluate the total flux {^/^^Iw) of both emission 
line components at each fiber from the quantities, / and w, presented in table 1 of Motta, 
Mediavilla, Munoz, Falco (2004). Error of flux at each fiber is estimated by using the error 
of the line intensity peak(57) and that of the line width {5w) as usual manner; the error is 
equal to ^/i^y/P ■ 6w'^ + • SP. Consequently, we can estimate the goodness of fit, i.e., the 
value, between a model and the observational data via equation 1. This procedure could 
be applicable if the emission line was detected, and if the decomposition into the broad line 
component and the narrow line component was succeeded. However, the emission line was 
not detected in some fibers, e.g., fiber 107 of Motta, Mediavilla, Munoz, Falco (2004). In 
such case, the total observed flux is set to be zero and the error is set to be the same as the 
error of a fiber that the total flux is the weakest in emission hne detected fibers. The value 
of such fibers can also be estimated via equation 1. Additionally, even if an emission line 
is detected, decomposition into the broad line component and the narrow line component 
might not be succeeded in some fibers, e.g., fiber 105 of Motta, Mediavilla, Munoz, Falco 
(2004). In such case, the emission line is treated as a single component in table 1 of Motta, 
Mediavilla, Munoz, Falco (2004), and the total observed flux is set to be an upper limit for 
both components and the error for both components is evaluated from the same observational 
data. If a model flux of a fiber exceeds the observed fiux of the fiber, the value of the 
fiber is evaluated by equation 1. If a model flux of a fiber is below the observed fiux of the 
fiber, the value of the fiber is set to be zero. 

The spatial sampling rate is set to be the same as INTEGRAL, and it is roughly 
0.5 arcsec in the case of Motta, Mediavilla, Munoz, Falco (2004). Both of the shape of 
a fiber cross section and the fiux loss between fibers were taken into account. As reported 
in Motta, Mediavilla, Munoz, Falco (2004), seeing is set to be 0.7 arcsec. We make use of 
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the data of 6 x 6 = 36 fibers, in total, in a rectangular array. Since the number of the 
fitting parameter is two; -Rnlr (or i^sLR for BLR) and A, degree of freedom in the fitting is 
36 — 2 = 34. The following results are summarized in the upper part of table 1. 

3.1. Broad line region 

The Gaussian intensity peak for the broad line component is presented in table 1 of 
Motta, Mediavilla, Munoz, Falco (2004), but the Gaussian width for the broad hne com- 
ponent is not. Thus, following the analysis by Motta, Mediavilla, Munoz, Falco (2004), we 
set the width to be 71 ± 11 A for all fiber. The searched parameter ranges of -Rblr and A 
are 50 ~ 2050 pc and 0.4 ~ 0.8, respectively. The best fit parameter search is performed 
with 50 pc resolution in -Rblr and 0.01 resolution in A. The result is presented in figure 3a. 
The best fit parameters are i?BLR = 50 pc and A = 0.47 (x^ = 28.01). However, even if we 
allow 1-a confidence level, a wide range of A values are acceptable and Rblr cannot be con- 
strained tightly. To get more reliable result, we had better to choose 3-a confidence region, 
rather than 1-a confidence region. Unfortunately, in such case, most of the parameter set is 
acceptable as seen in figure 3a, and we can say nothing about the lens model and the size 
of the BLR. Consequently, we cannot put any useful constraint on density profile of the lens 
model from the available data. We can only put an upper limit to the size of the BLR as 
-Rblr < 400 pc in 1-a confidence level. 

3.2. Narrow line region 

We use the same technique to investigate NLRs. The Gaussian intensity peak for the 
narrow line component and the Gaussian width for the narrow line component are taken 
from table 1 of Motta, Mediavilla, Munoz, Falco (2004). The searched parameter ranges and 
the resolution are the same as those in the case for the BLR. The fitting result is presented 
in figure 3b. Again, the best fit parameter is i?NLR = 50 pc and A = 0.47 (x^ = 25.09), 
but wider ranges of i?NLR and A are acceptable. However, there are two major differences 
between figure 3a and b. Firstly, the upper limit on /jInlr is significantly larger than that 
on it^BLR in 1-c confidence level. Secondly, the parameter ranges in the upper right part of 
figure 3b must be rejected, even if we allow 3-a confidence level. Unhke Motta, Mediavilla, 
Munoz, Falco (2004), we are not able to measure -Rnlr from the data. Its upper limit is 
~ 750 pc in 1-a confidence level. Again, we are not able to put useful constraints on the 
lens model as in the case of the BLRs. 
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4. Analyses by using Improved Data 

From the fitting results presented in the previous section, we understand that the cur- 
rently available data is not good enough to measure the size of the NLR and to constrain on 
the lens model. One reason for this can be attributed to the large spatial sampling rate in 
the used data, but this can easily be improved. Rather, better smaller spatial sampling rates 
with the same signal-to- noise (hereafter, S/N) ratio have already been achieved technically 
(Wayth, O'Dowd, Webster 2005). For instance, actually, GMOS-IFU on the Gemini North, 
SINFONI on VLT, or KY0T0-3DII on SUBARU telescope (Sugai et al. 2004) can provide 
us with two-dimensional spectra with 0.1 arcsec spatial sampling rate. Thus, it is worth of 
testing the potentiality of the proposed method for the data with 0.1 arcsec spatial sampling 
rate. In this section, we make mock observational data, which have better sampling rate, 
for a given parameter set (A and i?NLR or Rblr), and apply our method to those data to see 
how nicely it can reproduce the given parameter values. 

4.1. Our procedures 

Throughout this section, we adopt -Rnlr = 1000 pc and A = 0.5. In this lens model, 
the best fit value of c^e and e is 0.887 arcsec and 0.138, respectively. Then, we can calculate 
the expected observational images for a given spatial sampling rate, which is currently set to 
be 0.1 arcsec. Moreover, observational noise is artificially added to the expected image, and 
we obtain mock observational data. There may be several sources of noise, but we assume 
so-called photon noise to be dominant. Unless otherwise specified, seeing size is fixed to be 
0.7 arcsec ^ and the noise at the peak fiux is set to be 10 %. Data from 41 x 41 = 1681 
fibers are used for fitting. Roughly 4 arcsec x 4 arcsec field is covered, which is sufficient for 
our purpose. Since the number of parameters are 2 (i?NLR and A), degree of freedom for the 
fitting is 1681 — 2 = 1679. The following results are summarized in the lower part of table 1. 

4.2. Seeing effects 

The result is shown in figure 4a. The best fit parameter is -Rnlr = 1000 pc and A = 0.5 
(x^ = 1637.47). In 1-cr confidence level, the acceptable range of the parameters are still 

large, and the lens model degeneracy still remains in part. Even so, however, we can put 
tighter constraints on the parameters than for the currently available data (see figure 3) . The 
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acceptable ranges in 1-a confidence level are -Rnlr — 600 ~ 1300pc and A ~ 0.45 ~ 0.59. 
Therefore, we can nicely reproduce the given parameters with ~ 40% accuracy for i?NLR and 
~ 18% accuracy for A. 

The result in the case of smaller seeing size is shown in figure 4b. Again, the best 
fit parameter is -Rnlr = 1000 pc and A = 0.5 (x^ = 1634.92). Areas of all of the three 
confidence regions are smaller than those in the previous case (see figure 4a), although the 
differences are small. The acceptable ranges in l-a confidence level are -Rnlr — 650 ~ ISOOpc 
(~ 35% accuracy) and A ~ 0.45 ~ 0.58 (~ 16% accuracy). Although certain observations 
with better seeing will provide us with somewhat better information, the result do not alter 
much, compared with cases of the data with the same S/N ratio but with worse seeing. 

In contrast, if the incorrect seeing size is used in the analysis, the situation will become 
worse. The best-fit size is expected to be smaller (larger) than true value when we overes- 
timate (underestimate) the seeing size. For instance, if we assume 0.6 arcsec (0.8 arcsec) 
seeing and perform fitting to the mock observational data with 0.7 arcsec seeing, the min- 
imum value will exceed 3-a level, and the best-fit size and A will be 1500 pc (550 pc) 
and 0.46 (0.58), respectively. Thus, it must be crucial to know the seeing effect correctly 
for the success of our proposed method. The best-fit values are shifted along the sequence 
of degeneracy which has already been appeared in figure 4a. Since 0.1 arcsec corresponds 
to 850 pc at the redshift of this quasar, the difference between the true size and the best-fit 
size, ~ 500 pc, is somewhat smaller than the difference between the assumed seeing size and 
the true seeing size. This may be due to the most essential effect of gravitational lensing, 
image stretching of the observed images, in part. 



4.3. Observational errors 

Next, we examine the effects of observational errors on the results by decreasing noise, 
or increasing S/N ratio, at the peak ffux to 5%. In order to reduce observational errors to this 
value, roughly (0.1/0.05)^ = 4 times longer exposure time is needed than that in the previous 
case. The result is shown in figure 4c. Again, the best fit parameter is -Rnlr = 1000 pc 
and A = 0.5 (x^ = 1648.09), but each the confidence region becomes narrower dramatically. 
The acceptable parameter range in 1-a confidence level is -Rnlr — 850 ~ 1150pc (~ 15% 
accuracy) and A ~ 0.48 ~ 0.53 (~ 6% accuracy). Even in 3-a confidence level for this high 
S/N ratio case that is shown in figure 4c, the acceptable parameter ranges are somewhat 
smaller than those in 1-a confidence level for lower S/N ratio case that is shown in figure 4a. 
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4.4. Clumpy nature of the narrow line region 

In above studies, we assumed Gaussian emissivity profile for the NLRs as given by 
equation 2. In actual Seyfert galaxies, however, the emissivity profile may not be smooth 
but can be very complex; likely to be composed of small numerous clumps. Here, we assume 
that the NLRs are composed of 100 clumps, the emissivity profile of each clump obeys 
Gaussian profile (equation 2) with a size oi R — 50 pc. We assume that the clumps are 
randomly distributed and the distribution of distance between the source center and each 
clump obeys Gaussian distribution with the typical size of 1000 pc. 

The results are shown in figure 4d. We find that the best fit values of the parameters 
are i?NLR = 1000 pc and A = 0.49 (x^ = 1660.11). The best fit value of A is slightly smaller 
than the given value, but that of i?NLR is identical to the given value. The shape of 1-a 
confidence region is different from that in figure 4a, but the overall properties, such as the 
direction of the major axis of 1-a confidence region, are similar to each other. Although 
the extent of 1-a confidence region in this figure is smaller than that in figure 4a, the given 
values are included in 1-a confidence region, i?NLR = 650 siml200 pc and A = 0.45 ~ 0.56. 

Therefore, we can safely conclude that we can put a reliable constraint on the size of the 
NLRs and the lens model by the proposed method, even if the clumpy nature of the NLRs 
are taken into account. The clumpy nature may work as an origin of the systematic error in 
the proposed method. 

5. Concluding Remcirks 

In this paper, we present the potentially useful method to constrain the size of the 
NLR and the lens model. Taking the case of Q2237+0305 as an example, we re-analyze 
the two dimensional spectroscopic data taken by Mediavilla et al. (1998) and obtain reliable 
constraints on the size of the BLR and the NLR. Although both of the sizes are consistent 
with point-like source, the possibility of the extended source is still remains. The obtained 
upper hmit to the sizes of the BLR and the NLR are 400 pc and 750 pc, respectively. As 
far as the mass profile of the lens galaxy is nicely represented by a model where the lens 
potential is elliptical and a power-law in radial profile, as we assumed in this paper, these 
values can be appropriate upper limits. 

In addition, we have demonstrated that the currently available instruments can provide 
a better opportunity to explore the NLRs and the lens model. Seeing is an important 
factor in ground-based observations. However, it may not seriously affect our results, if we 
know the seeing effect correctly. Observations with higher S/N ratios and with higher spatial 
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sampling rates will be useful to give better constraints on the parameters. Even if the spatial 
distribution of the NLR is not smooth but clumpy, we will be able to reproduce the values 
of the basic parameters successfully. 

One thing that we have to mention is so-called "mass sheet degeneracy" in gravitational 
lensing phenomena; if there is a mass sheet with the surface mass density of Ks between 
observer and the source, any observed structure will be scaled up by a factor of (1 — Ks)~^ 
in length and by a factor of (1 — Ks)~^ in area. Consequently, the actual source size should 
be the produce of (1 — /t) and the obtained source size ^. Though the measurement of Kg 
is practically hard and annoying issue, such difficulty may be solved at least in part by 
observational studies of the environments around lens systems (e.g., Faurc ct al. 2004). 

If the emissivity profile of the NLR has an elliptical shape or a double-cone shape, the 
lensed image of the NLR can be different from those considered in this paper and these 
different properties will alter the results to some extent. Further, the existence of substruc- 
tures in or around the lens galaxy may also affect the observed image of the NLR (e.g., 
Metcalf, Moustakas, Bunker, Parry 2004). To evaluate such effects, quantitative estimations 
are necessary. Even in such case, however, we can, in principle, find correct values of the 
parameters in a statistical fashion, e.g., by searching the deepest valley of (see figure 4), 
and to obtain an observational insight into the shape of the NLR and/or substructures as 
well as the size of the NLR and the lens model. 

The author acknowledges to S. Mincsliigc, H. Sugai, T. Nagao, K. Yahata, W.-H. Bian, 
and anonymous referee for their valuable comments and discussions. 
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Table 1. Summary of the fitting results 



data 


size^ 


(1-(T range) 


A 


(l-(j range) 


figure 


BLR^ 


50 


(50 ~ 400) 


0.47 


(0.42 ~ 0.71) 


figure 3a 


NLR^ 


50 


(50 ~ 750) 


0.47 


(0.42 ~ 0.79) 


figure 3b 


standard^ 


1000 


(600 ~ 1300) 


0.50 


(0.45 ~ 0.59) 


figure 4a 


better seeing'^ 


1000 


(650 ~ 1300) 


0.50 


(0.45 ~ 0.58) 


figure 4b 


smaller errors'^ 


1000 


(850 ~ 1150) 


0.50 


(0.48 ~ 0.53) 


figure 4c 


clumpy nature'^ 


1000 


(650 ~ 1200) 


0.49 


(0.45 ~ 0.56) 


figure 4d 



^-Rblr for the BLR, and -Rnlr for the NLR (unit: pc). 

^Real observational data obtained by Motta, Mediavilla, Munoz, Falco 
(2004) (sec section 3). 

■^Mock observational data calculated in this paper (see section 4). 
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Fig. 1. — Total value (the bottom panel), the best fit value of a-^ in the unit of arcsec 
(the middle panel), and that of e (the top panel) is presented with solid line as a function A. 
In the bottom panel, 1-, 2-, and 3-cr confidence level is also presented by the dotted, dashed, 
and long-dashed line, respectively. In the top panel, a critical value of e which is noted by 
Kassiola, Kovner (1993) (see section 2.2) is also presented with the dashed line. 
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Fig. 2. — Calculated images that are to be compared with observational data. We indicate 
the values of -Rnlr and A in the bottom and in the left of the panels, respectively. The flux 
of each pixel is normalized by the peak flux. We assume 0.7 arcsec seeing and 0.1 arcsec 
sampling rate. 
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Fig. 3. — The results of fitting to the broad hne component (a: left panel) and the narrow 
line component (b: right panel), respectively. The thick, middle, and thin lines present 1-, 
2-, and 3-cr confidence levels, respectively. The filled circle indicates the best fit parameter. 
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Fig. 4. — Same as figure 3, but for the moclc observational data (see section 4). Panel (b) 
and (c) shows the result for better seeing case and smaller error case compared with panel 
(a), respectively. Panel (d) shows the result for the clumpy NLR. 



